Metamaterials proposed as perfect magnetoelectrics 



A. M. Shuvaev,^ S. Engelbrecht/ M. Wunderlich,^ A. Schneider/'^ and A. Pimenov^ 

^ Experimentelle Physik IV, Universitdt Wiirzburg, 97074 Wilrzhurg, Germany 
^Fachbereich 1, Universitdt Bremen, 28359 Bremen, Germany 
(Dated: April 27, 2010) 

Magnetoelectric susceptibility of a metamaterial built from split ring resonators have been inves- 
tigated both experimentally and within an equivalent circuit model. The absolute values have been 
shown to exceed by two orders of magnitude that of classical magnetoelectric materials. The meta- 
material investigated reaches the theoretically predicted value of the magnetoelectric susceptibility 
which is equal to the geometric average of the electric and magnetic susceptibilities. 



Magnetoelectric effect manifests a connection between 
electricity and magnetism and has been predicted already 
1860 by Pierre Curie, but could be confirmed experimen- 
tally only about half a century ago [TJ [2]. Up to now 
the experimentally observed effects remain rather weak. 
The weakness of the magnetoelectric susceptibility can to 
some extent be explained by the absence of strong mech- 
anisms which couple magnetism and electricity on the 
microscopic level. As has been shown theoretically [HIS], 
the allowed limiting value of the magnetoelectric suscep- 
tibility (^) is quite large and equals the geometric average 
of electric (xe) and magnetic (xe) susceptibilities: 

< XeXm (1) 

In classical magnetoelectric materials like Cr2 03 the lim- 
iting value of Eq. ([T]) is failed by about two orders of 
magnitude ^ ,3^ . In efforts to increase the value of the 
magnetoelectric effect, materials revealing both strong 
electric and magnetic susceptibilities have been brought 
into consideration. Especially close to phase transitions, 
the electric and magnetic susceptibilities may diverge 
in ferroelectrics and ferromagnets. Materials simultane- 
ously showing the ferroelectricity and ferromagnetism are 
called multiferroics and they are presently the subject of 
intensive research |2l |41 . 

Assuming that the equality in Eq. ([T]) holds, it can be 
rewritten in the form ^ = s/XeXm and the constitutive 
relationships 

M = XmH^XmeE (2) 
P = XeE^XemH , (3) 

can be reduced to 

M = iP,f^ , (4) 

V Xe 

i.e. electric polarization and magnetic moments must be 
directly proportional to each other. Here we have rewrit- 
ten Xem = and utilized the symmetry of the magneto- 
electric coefficients Xme = —Xem niE]- The condition in 
Eq. Q has been assumed in the first molecular theories 
for magnetoelectric effect [1 nearly a century ago. 



We note that within the arguments of thermodynamic 
stability [6 another magnetoelectric inequality can be de- 
rived: < ejii which utilizes the magnetic permeability 
/i = 1 + Xm and electric permittivity e = 1 -\- Xe instead 
of susceptibilities in Eq. ([T]). 

Due to their unique electrodynamic properties, meta- 
materials (i.e. artificial materials) [7-9 may show new 
ways to solve the problem of the weakness of magneto- 
electric coupling. Some examples of breakthroughs in 
various topics of modern electromagnetism which were 
stimulated by metamaterials are: reversing the laws of 
conventional optics [10 , cloaking [llill2j, or overcoming 
the resolution limit of optical devices [I3l [14] . Mixing of 
electric and magnetic responses is another useful prop- 
erty of metamaterials which can be utilized to generate 
new effects. As has been shown recently, the magneto- 
electric coupling in metamaterials lead to strong optical 
activity p!5HT9] which is directly connected to intrinsic 
chirality. In agreement with the theoretical prediction 
[2Q[ I2TI the chirality of metamaterials can lead to giant 
polarization rotation and provide another routes to ob- 
tain negative refraction [22-:23^ for circularly polarized 
waves. Recently, for metamaterials with zero permittiv- 
ity and permeability an estimate for the limiting value 
of the magnetoelectric effect Re{^) < Im^/eJI has been 
obtained [21]. 

Many designs of metamaterials are based on split ring 
resonators [26 . These elements can be seen as the small- 
est possible representations of the well known LC-circuit 
with a single inductance loop as given by the metallic 
ring and a tiny capacitance produced by the gap in the 
ring [27 . Split ring resonators have been originally de- 
veloped to achieve negative magnetic permeability which 
is a key property for design of metamaterials with neg- 
ative refraction [7 . As has been realized recently [28j, 
the split ring resonators strongly modify the interac- 
tions with electromagnetic radiation by introducing a so 
called bianisotropy term into the set of basic equations 
[29] , which is closely connected to the magnetoelectric ef- 
fect [2j [5] . This additional term in the constitutive rela- 
tions cross-couples the magnetic and electric fields within 
a split ring resonator [30l |3]j . The bianisotropy offers an- 
other degree of freedom [32' in controlling the properties 
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of light. Here we note that, contrary to chiral metama- 
terials, the bianisotropy does not automat icahy lead to 
polarization rotation for geometries parallel to the princi- 
pal optical axes. In order to obtain polarization rotation, 
these structures must be tilt or measured within off-axis 
geometry. The corresponding effects [33l [34] have been 
termed extrinsic chirality. 

In this work we show that metamaterials built from 
split ring resonators achieve magnetoelectric effects equal 
to the theoretically limiting value in Eq. ([T]). To prove 
this we investigate a metamaterial of split ring resonators 
within different geometries, especially including those 
sensitive to magnetoelectricity. This allowed to obtain 
electric, magnetic and magnetoelectric susceptibilities 
and compare them with a simple circuit model. We show 
that the metamaterial investigated indeed reaches the 
theoretical limit for magnetoelectric coupling. This value 
is due to direct proportionality of electric and magnetic 
moments in split ring resonators. 

Transmittance experiments at millimeter-wave fre- 
quencies (60 GHz < u < 120 GHz) were carried out in a 
Mach-Zehnder interferometer arrangement [35]|36]. This 
arrangement allows to measure both the intensity and 
the phase shift of the radiation transmitted through the 
sample within controlled polarization geometries. The 
split ring resonator arrays used in present experiments 
were prepared by chemical etching of copper-laminated 
board. The rings are typically 0.35 mm x 0.35 mm in size 
with the gap width d = 0.17 mm. The lattice constant of 
the metamaterial is / = 0.7 mm. The characteristic pa- 
rameters of various sets of the split ring resonators have 
been varied within approximately a factor of two and 
showed qualitatively similar results. Woven glass with 
a thickness of 0.56 mm was used as a non-conductive 
substrate. The refractive index of the substrate has been 
determined in a separate experiment as Ug = 2.07+0.04z. 

Split ring resonators seem to represent an ideal mag- 
netoelectric material fulfilling the condition XmeXem = 
XeXm- Indeed, from the effective RLC-circuit model and 
simple calculations we get 

Xe = nC-d^- F{uj) Xm = nC-S^'^- F{uj) (5) 



Xe 



-Xme nC • dS^ • F{uj) 
c 



(6) 



Here n is the density of the rings, d and C are the effective 
width and capacitance of the gap, S is the area of the 
rings, and uo is the angular frequency. We use the Lorentz 



substitution F{uj) = uJq/{uJq 



iujj) where ujo and 



7 are resonance frequency and width, respectively. The 
symmetry of the magnetoelectric coefficients is fulfilled 
automatically in this model. 

In order to obtain the electrodynamic parameters of 
the split ring resonators we have carried out the ini- 
tial experiments within the experimental geometries sug- 
gested in Ref. [29 (shown in the insets to Fig. [T]). 
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FIG. 1: Transmission characteristics of the metamaterial of 
split ring resonators. Two experimental geometries shown 
are mostly sensitive to the dielectric (green) and magnetic 
(red) contributions. Symbols - experiment, lines are fits us- 
ing Lorentzian characteristics of the split rings. Solid lines - 
fits include magnetoelectric coupling, dashed lines - magneto- 
electric susceptibility is set to zero, demonstrating negligible 
influence on the spectra. The pictograms show the geometry 
of the experiments. 



Three relevant geometries in this case are "magnetic" 
{h perpendicular to the plane of the rings), "electric" 
(e parallel to the gap of the rings) and "magnetoelec- 
tric" (both excitations are realized simultaneously). In 
these three geometries and within reasonable approxi- 
mation the effective refractive indexes basically deter- 
mine the transmittance close to the resonance and they 
are given by: Ue = V^^^^M^ = - ^^/s:, and 
^me = — respectively. Here we neglect the influ- 
ence of the substrate for simplicity. Although the mag- 
netoelectric susceptibility is included in these equations, 
the dominating terms for typical parameters of the model 
are given by a/^, and ^/eJI. In all cases the mag- 
netoelectric susceptibility represents a weaker correction 
under the square root. This is demonstrated in Fig. [l]as 
we set the magnetoelectric susceptibility to zero, which 
simply leads to a shift of the resonance frequency. In all 
series of experiments with varying geometries the influ- 
ence of the magnetoelectric susceptibility was below the 
experimental accuracy. This accuracy depends not only 
on experimental uncertainties, but also on the assump- 
tions of the circuit model, like neglecting of the cross cou- 
pling effects, or assumption of infinitely small sizes of the 
rings compared to the wavelength. On the contrary, elec- 
tric and magnetic geometries robustly depends on electric 
(Xe) and magnetic (Xm) susceptibilities. Therefore, both 
susceptibilities may be determined from the spectra in 
Fig.B 

The result of the experiments described above may now 
be extended to obtain the magnetoelectric susceptibility 
using further geometries of the experiment. In order to 
get better sensitivity to the magnetoelectric susceptibil- 
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FIG. 2: Transmission and phase shift spectra of the metama- 
terial within various tilted geometries as indicated. Symbols - 
experimental data within the geometries as indicated, circles 
correspond to the spectra with parallel polarizers, triangles - 
with crossed polarizers. Solid lines are fits within the 4x4 
matrix formalism 38 and using Eq. (5[6) for susceptibili- 
ties curves. In the panel C the model curves with crossed 
polarizers coincide for both geometries. 



ity, the sample must be measured in the tilt geometry 
and the signals in parallel and in crossed polarizers have 
to be compared (see inset in Fig. [2| . An example of such 
tilt geometries has been presented in Ref. ^31j , where the 
existence of a cross-coupling terms has been detected. 
The effectiveness of the tilt experiments can be demon- 
strated within the simplified assumptions. If we assume 
ka << 1, then the corresponding Maxwell equations may 
be solved more easily leading to analytical expressions for 
all relevant geometries (e.g. Ref. [37]). Here /c = 27r/A is 
the wavevector of the electromagnetic wave and a is the 
sample thickness. If we further neglect the (/ca)^ terms 
compared to the terms linear in /ca, then e.g. for both 
geometries in the left panels in Fig. |2] and within 45° 
incidence the following expression for the complex trans- 
mittance in crossed polarizers can be written: 



tito) 



ka • ^ 



2V2/i + ika[l + ^2 _ ^ 2)] 



(7) 



The last formula clearly shows that the transmittance in 
crossed polarizers is directly proportional to the magne- 
toelectric susceptibility. Therefore, this geometry pro- 
vides a sensitive tool for magnetoelectric effect. In the 
following, the calculations have been done within the ex- 
act 4x4 matrix formalism as described in Ref. [38]. 

Figure [2] shows transmittance and phase shift spectra 
in four most important tilt geometries, both in parallel 
and in crossed polarizers. The main result of this data 
is given in the left bottom panel of Fig. |2] showing the 
nonzero transmittance in crossed polarizers close to the 
resonance of the rings (109 GHz). As expected already 
from the simplified equation Eq. ([7|), reasonably strong 
signal can be observed in geometries with crossed polar- 
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FIG. 3: Permittivity, permeability and bianisotropy of the 
metamaterial used in the present work. Symbols - experiment, 
lines - circuit model given in Eqs. ([5|[6]). The magnetoelectric 
susceptibility has been assumed to obey ^1 ^XeXm — 0.9 =b 
0.15 instead of Eq. 



izers (red and green triangles). In the left bottom panel 
of Fig. |2] the nonzero signal below about 90 GHz for 
crossed geometry are due to the radiation leakage around 
the sample and should be neglected. The leaky signal 
can be also seen in the phase shifts data as shown in the 
left upper panel. Below about 90 GHz the amplitude of 
the signal is getting low and the phase measuring system 
looses the stability (red and green triangles). 

We note another interesting feature of Fig. |2j In the 
geometry with ac electric fields within the plane of the 
rings and parallel to the gap (black circles in the right 
panels of Fig. [2| the resonance is clearly seen in paral- 
lel polarizers. However, we observe no signal in crossed 
polarizers in this geometry (black triangles). In two ge- 
ometries presented in the right panels of Fig. |2] solely a 
stray signal could be detected. This result is supported 
both by simple and rigorous theories. Both calculation 
methods predict zero amplitude of the perpendicular po- 
larization and reflect the symmetry of the problem. 

The results of the transmittance experiments in per- 
pendicular (Fig. [1]) and tilted (Fig. |2| orientations are 
sufficient to determine all three complex susceptibilities 
for the metamaterial of split ring resonators. The re- 
sults of the self-consistent calculation of these parame- 
ters based on the 4x4 matrix formalism is presented in 
Fig. [3] As could be already expected from the fit of the 
transmittance spectra, all three susceptibilities reveal a 
resonance-like form. The results in Fig. |3] clearly demon- 
strate that the metamaterial of split ring resonators in- 
deed models a perfect magnetoelectric with = XeXm- 
Here solely a factor of 0.9 ± 0.15 had to be introduced 
in order to obtain the self-consistence of the data. Tak- 
ing into account the assumptions made, this factor most 
probably reflects the uncertainties of the experiment. Fi- 
nally, we recall that the condition = XeXm is expected 
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for all materials revealing direct proportionality between 
electric and magnetic moments (M = aP). Here a is a 
material constant which equals a = iSuj/cd for split ring 
resonators. 

In summary, using millimeter wave spectroscopy of 
the complex transmission coefficient the electrodynamic 
properties of a metamaterial of split ring resonators have 
been investigated. The sensitivity to the magnetoelec- 
tric effect has been obtained within tilt sample geome- 
try and calculated within 4x4 matrix formalism. We 
prove experimentally and within a circuit model calcula- 
tion that metamaterials from split ring resonators reach 
the maximum theoretical values of the magnetoelectric 
susceptibility limited by < XeXm- This value appears 
to be about two orders of magnitude above the typical 
coupling constants for conventional magnetoelectrics like 
CrsOs. 
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